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Abstract 

This current work is an extension to work previously done by 
the authors. The dynamics of a baby skyrmion configuration, in a 
model Landau-Lifshitz equation, was studied in the presence of vari- 
ous potential obstructions. The baby skyrmion configuration was con- 
structed from two Q = 1 hedgehog solutions to the new baby Skyrme 
model in (2 + 1) dimensions. The potential obstructions were created 
by introducing a new term into the Lagrangian which resulted in a 
localised inhomogeneity in the potential term's coefficient. In the bar- 
rier system the normal circular path was deformed as the skyrmions 
traversed the barrier, after which the skyrmions orbited the boundary 
of the system. For critical values of the barrier height and width the 
skyrmions were no longer bound although the unbound behaviour is 
not clearly distinct from the bound. In the case of a potential hole the 
dynamics of baby skyrmions is dependent upon the binding energy 
of the system. Depending upon its value, the skyrmions' behaviour 
varies. The angular momentum must be modified to ensure overall 
conservation. We show that there exists a link between the oscillation 
in the skyrmion's energy density and the periods of non-conservation 
of the angular momentum in Landau-Lifshitz models. 
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I Introduction 



There has been recent and considerable interest in the scattering of topo- 
logical solitons off potential obstructions in a variety of different models. 
Al-Awali et al pQ, [2], [3] have considered the scattering of various topolog- 
ical and non-topological solitons off potential holes and barriers. Speight 
[TT] has recently investigated the scattering of a topological soliton in a anti- 
ferromagnetic system, where the continuum dynamical equation resembles a 
second order relativistic wave equation. Both of these system have shown 
that the scattering of topological soliton off holes or barriers exhibit some 
interesting results. Topological solitons are, of course, classical objects. How- 
ever, they describe extended objects and, as shown in [I], some of their prop- 
erties resemble those of quantum systems. We have recently investigated the 
scattering properties of a baby skyrmion configuration in a magnetic system 
[5], whose motion is then governed by the Landau-Lifshitz equation. This 
system showed some interesting and non-intuitive results. The potential ob- 
structions were placed symmetrically about the x-axis for all values of x. 
In this current study we shall investigate the same system with a different 
obstruction geometry. 

The baby skyrmions we are interested in are those of the (2 + 1) new baby 
Skyrme model. The Lagrangian for the (2 + 1) new baby Skyrme model is 
given by: 

c = jndp£ ■ - \ l2 [{d^ ■ d^f - (d^ ■ dA){d»± ■ - ~ 73 (i - 4), (i) 

where 71,2,3 are positive constants and is a triplet of scalar fields: <p(x, t) = 
{<pi(x, t);i = 1, 2, 3}. The fields of this model are subject to the constraint 0- 
0=1. The latter two terms have been introduced to avoid the consequences 
of Derrick's theorem [6] and to stabilise topological soliton solutions in two 
dimensions. The condition that 2 = 1 is imposed so that the target space 
is the 2-sphere, such that is now a map : M? — » S 2 . For finite energy 
solutions it is necessary for the fields to tend to a vacuum at infinity, where 
03 = 1 at 00. This results in a compactification of M 2 so that now takes 
values in the extended plane M 2 U 00, which is topologically equivalent to 
S 2 . The constraint equation 2 = 1 and the boundary condition at infinity 
results in the field becoming a non-trivial map : S 2 -> S 2 . Each soliton 
solution is grouped into a different homotopy class according to the winding 
number, or topological charge, of this map. The topological charge Q is given 
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by: 

Q = ±- [ d 2 xe ij( p- (drfxdrf), (2) 

OTT J R 2 - - - 

where the indices i,j run over the space coordinates and Q G Z. The 
topological soliton solutions of the new baby Skyrme model are called baby 
skyrmions. Baby skyrmion solutions for this model are constructed using 
the Hedgehog anzatz and must be found numerically. This problem reduces 
to solving a second order differential equation for the profile function f(r). 
The configuration we use in our investigations consists of two single charged, 
Q = 1, baby skyrmions. The construction of these solutions was described 
in our previous work [5]. Static solutions of the new baby Skyrme model are 
static solutions of the Landau-Lifshitz model. It has been shown that the 
dynamics of two interacting baby skyrmion solutions in a Landau-Lifshitz 
model, resemble the behaviour of the experimentally observable magnetic 
bubbles [9], [10]. The Landau-Lifshitz equation and its constraint are given 
by: 

i = *x- ¥ , f=l, (3) 
where E is the energy functional written as: 

E= ! d 2 x8, 



and S is the static energy density of the new baby Skyrme model given by: 

£ = \-ndi± ■ 04 + ■ d^f - {84 ■ d$(d4 ■ d&] + ^73(1 - 

Analysis of the dynamics in Landau-Lifshitz systems has been greatly 
simplified by the work of Papanicolaou and Tomaras [7], who constructed 
unambiguous conservation laws for the system governed by ([3]) . In their work 
they found that the important quantity was the topological charge density 
q: 

q = e ij $- (dj-0x (4) 
Some of the conservation laws can be constructed as a moment of q. They 
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involve: 

I — - f d 2 xx 2 q, (5) 

m = [ d 2 x(<p 3 -l), (6) 
J = l + m, (7) 

where / is the orbital angular momentum, m is the total magnetization in 
the third direction and J is the total angular momentum. Conservation laws 
for the system were constructed by examining the time evolution of q: 

q = -eijdididjt, (8) 

where d\Oji can be written in terms of the energy functional E: 

fail = ■ 9 A ■ ( 9 ) 



Taking an explicit time derivative of ([5]) gives: 

i = - / d 2 xx 2 q. (10) 
2 Jk 2 

The guiding centre coordinate R of the soliton is defined as the first moment 
of the topological charge density q: 



R=-— d'xxq. (11) 
4vrQ J R 2 

In our construction of the potential obstructions we adopt a similar ap- 
proach to [1] , [5] and introduce a term into the Lagrangian ([T|) which vanishes 
in the vacuum state 03 = +1. The obstructions need to be introduced in this 
manner so the tails of the skyrmions are not changed by the obstruction. 
Therefore, as in [5], the additional potential term added to the new baby 
Skyrme Lagrangian (JTJ), modifies the potential coefficient such that it now 
depends on the space coordinates. The introduction of this term implies that 
the static part of ([T|) can be rewritten as: 

£ = ipidd-d^- ^72 \{d4- dif) 2 - (drf-djcjj) (<%!>• djcjj)] - -j 3 (x,y)(l-(f>l), 

(12) 
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where the potential term coefficient 73 is now a function of the coordinates 
(x, y) and the static part of (pQ) is only considered as imposed by the Landau- 
Lifshitz equation. In our previous work we presented the results of our inves- 
tigations into the scattering properties of a baby skyrmion configuration off 
various symmetric potential obstructions. We can ask ourselves, how many 
of the properties of that system were indicative of its symmetric nature? In 
other words, if the symmetric nature of the system was to disappear, will the 
dynamics of the baby skyrmions be greatly changed? Hence in this paper we 
shall discuss the results and analysis of the same baby skyrmion configuration 
used in [5], in the presence of an asymmetric potential obstruction of width 
b and height or depth T. The asymmetric obstruction is introduced into the 
system in the same manner as the symmetric obstruction where 73(2;,?/) of 
Q12p in this system can be written as: 

73 (x, y) = i + r{e(-x + x ) [Q(y + yo) - <s>(y - y )] }, 

where yo = b/2 and xq = 0. The sign of the additional potential term's 
coefficient, T, determines whether the obstruction is a hole, T < 0, or a 
barrier, T > 0. 

The evolution of the dynamics of the baby skyrmion configuration, in- 
teracting with the potential obstructions, has to be done using numerical 
methods. The details of the numerical procedures used in integrating the 
Landau-Lifshitz equation or the discretisation of the continuum problem etc 
are discussed in the appendix. In the absence of any obstructions, T = 0, 
the two baby skyrmions of the configuration orbit each other along a circular 
trajectory. The baby skyrmions orbit the centre of the configuration (0, 0) 
anti-clockwise. In the rest of the paper we shall refer to baby skyrmions as 
simply skyrmions. 

II Potential barrier 

In this section we shall illustrate some of the results seen in the potential 
barrier system for various values of the obstruction height T and the barrier 
width b. The physically interesting part of these systems is the initial phase of 
the dynamics, when the skyrmions are alternating between moving 'on' and 
'off' the barrier. The motion here is very similar to that seen in the symmetric 
system of [5] . Due to the asymmetry of the system the skyrmions no longer 
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Table 1: Table showing the variation in the total energy and binding energy 
for a d = 6 two skyrmion configuration interacting with a potential barrier 
of width b = 2 for various values of T . 



execute a symmetric path about the obstruction. Once the skyrmions have 
drifted away from the potential barrier they move along the boundary of the 
system. 

The analysis of the skyrmion scattering in the barrier system can be 
divided into the sections where the skyrmions are far from the obstruction 
and when they are executing their circular motion 'on' or near the barrier. 
The initial phase of the dynamics corresponds to the latter of these two 
sections. Fig. [T] shows multiple plots of a d = 6 skyrmion configuration 
interacting with an asymmetric barrier of width b = 2 and depth T = 0.25. 
The skyrmions centres are initially located at (0, ±3) where we refer to the 
skyrmion initially placed at (0, 3) as the upper skyrmion and (0, —3) as the 
lower one. Fig. dh and b show the trajectories of the upper and lower 
skyrmions over a time length of 3500 sees. Fig. [Tb shows the two trajectories 
plotted together. It can be seen from the first two plots that initially the two 
skyrmions try to execute their normal circular motion around each other. 
As the skyrmions approach the obstruction, their path is deformed due to 
the barrier. The trajectory of the skyrmions in this system suggests that 
the skyrmions still form a bound state for these values of T and d. The 
binding energy of a two-skyrmion configuration in the presence of a potential 
obstruction is denoted by Eb and can be defined by the following: if the 
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a b 

Figure 1: Plots of the trajectories of the upper (a) and lower (b) skyrmions 
in a system with an asymmetric potential barrier of width b = 2 and height 
T = 0.25, where the horizontal axis represents the x coordinate and the 
vertical axis the y coordinate. 
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energy of the skyrmion configuration in the presence of a barrier is denoted 
by E 2 and E\ is the energy of a single skyrmion placed at the same position 
as one of the skyrmions in E 2 , then the binding, or the interaction energy, is 
given by: 

E B = E 2 - 2E X . (13) 

Table [T] shows the variation in the binding energies for various values of V 
for a d — 6 skyrmion configuration, interacting with a potential barrier of 
width 6 = 2. Here for T = 0.25 the binding energy of the skyrmions Eb < 
and thus the skyrmions still form a bound state, which the plots in Fig. [1] 
confirm. 

Fig. [2] shows the plot of the initial phase of the above system, up to 1500 
sees, with a set of labels A-K and 1-3. The labels A-K are reference points 
used to examine the motion and speed of the skyrmions during this initial 
phase. The labels 1-3 refer to the different times the skyrmions traverse the 
barrier. The potential obstruction's position is outlined by the dotted lines 
appearing on the plot of Fig. [2j Although the barrier has been defined in 
section fl] the lines have been continued for x > and for y > y , y < —y 
for purely illustrative purposes. In order to identify the times at which the 
skyrmions have entered the barrier region we have further plotted x(t) and 
y(t) for the upper skyrmion shown in Fig. [31 These plots are labelled with 
A-K accordingly. 

The initial path of the upper skyrmion as it approaches the barrier is 
A-B. The distance of A-B and the distance of the path over the barrier B-C 
are comparable. However, the time taken for the skyrmion to travel A-B is ~ 
100 sees, whereas to traverse B-C it is ~ 30 sees. If we examine the gradient 
of y(t) for the upper skyrmion, we can see that the gradient rapidly increases 
and then decreases as the skyrmions traverse the barrier. The speeding up of 
the skyrmions as they traverse a potential barrier was seen in the analysis of 
systems with a symmetric potential obstruction. In the further two scattering 
paths the times taken to traverse the barrier is ~ 30 sees for F-G, and ~ 30 
sees for J-K. It can be shown that the lower skyrmion traverses the barrier 
in similar times as the upper skyrmion, for paths 1, 2 and 3. 

Fig. H] shows a plot of the skyrmion separation as a function of time, d(t), 
for the first 1500 sees, with an initial separation of d — 6. As the skyrmions 
move 'on' and 'off' the barrier, the distance of separation decreases and 
increases accordingly. The minima of d(t) correspond to the times when 
one of the skyrmions is on the barrier. The skyrmions, initially placed in 
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Figure 2: Plot of the trajectory for the upper skyrmion of a d — 6 con- 
figuration interacting with an asymmetric potential barrier of width b = 2 
and T = 0.25, where the horizontal axis represents the x coordinate and the 
vertical axis the y coordinate. 



the attractive channel, compensate for the increase in potential energy by 
reducing their distance of separation. Once off the barrier the skyrmion 
configuration returns to its initial state and the skyrmions separate once 
again. This process repeats with each of the different paths 1-3. 

II. 1 Away from the barrier 

Another interesting feature of this system is the way the skyrmions avoid the 
boundary and potential obstruction after the initial phase described previ- 
ously. After the transient period where the skyrmions traverse the barrier, 
they eventually settle down to a steady state where the skyrmion configura- 
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Figure 4: Plots of the distance of separation, d, between two skyrmions in a 
system with an asymmetric potential barrier of width b = 2 and T = 0.25. 



tion moves about the grid avoiding the barrier and the boundary. During this 
steady state the skyrmions continue to execute the circular motion about the 
configuration centre when they are close enough to interact. The skyrmions 
have a strong repulsion from the boundary and keep a certain distance from 
it. The skyrmions' motion in this phase is determined by the boundary. If 
we increase the size of the boundary, after performing a similar initial phase, 
the skyrmions continue to move around the edge of the boundary. This is 
illustrated in Fig. [5] which shows the upper skyrmion of the same system of 
b = 2 and T = 0.25, plotted with an identical upper skyrmion in a system 
with the same values of b and T but with a grid size of 291 x 291 points. 
The boundary of the two different grid sizes is also shown on Fig. [51 The 
larger grid means that the boundary of the system is further away than for 
the system with a smaller grid size. This plot confirms that the skyrmions' 
trajectory away from the barrier is dependent upon the boundary. Thus 
the physically relevant part of the dynamics i.e. the phase unique to the 
problem of interest, is the initial phase before the skyrmions drift along the 
boundary's edge. 
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Figure 5: Plot of the trajectories of the upper skyrmions, where both 
skyrmions are initially at (0,3), in a system with an asymmetric potential 
barrier of width b = 2 and T = 0.25 for grid sizes 251 x251 and 291 x 291 
with the boundary of different grid size shown. The horizontal axis represents 
the ^-coordinate and the vertical axis the y-coordinate. 
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In the previous sections we described the general system of a bound 
skyrmion configuration interacting with a potential barrier. This system was 
chosen because it illustrated many of the generic features of the skyrmion's 
motion in these systems either in the region near the obstruction or away 
from it. However, there are a few additional features which exist for only 
certain values of the barrier height T, and distance of skyrmion separation 
d, and shall be presented here. 

In the d = 6, T = 0.25 barrier system, we observed that the skyrmions 
were still bound due to value of the binding energy. This was also reflected in 
the trajectories of the skyrmions in that system. The skyrmion configurations 
when d < 6, due to the attractive channel, are more tightly bound than the 
d = 6 skyrmions. In systems with small values of T, the skyrmions exhibit 
bound behaviour in both the binding energy and trajectories. Away from 
the barrier they orbit the boundary of the system as usual. This behaviour 
is evident in Fig. [6] which shows the plots of the trajectories for a skyrmion 
configuration with d = 4 and d = 5 for various values of the barrier height 
T. The skyrmions in these systems orbit each other much faster than the 
looser bound d = 6 system. The trajectories of the skyrmions in the d = 4, 5 
systems are slowly distorted as they execute their motion. As the value of T 
increases the distortion of the skyrmion trajectory becomes more pronounced. 

In the absence of any obstructions the skyrmions form a bound state. 
The strength of the bound state is dependent upon the separation between 
the two skyrmions. For d = 6 the skyrmions are reasonably well separated 
and the skyrmions are not tightly bound. For large values of T however, 
the binding energy of these skyrmion configurations can increase and the 
bound state can be broken. As the barrier height gets larger, the skyrmions 
are unable to overcome it. In [5], in the barrier's transition dynamics, it 
was seen that the skyrmions cannot distinguish between the boundary of the 
system and very large barriers. We showed that the skyrmions scatter at an 
angle a relative to the barrier's edge. This angle a — > as V gets larger. In 
the asymmetric system for a sufficiently large barrier, which the skyrmions 
cannot traverse, the skyrmions begin to move at this angle a. Initially the 
upper skyrmion starts to move in the direction away from the obstruction 
but due to the asymmetry of the system it is able to interact with the lower 
skyrmion and sling shot around the edge of the barrier. 

The trajectory of a set of skyrmions, where the bound state is broken by 
the obstruction, is much more difficult to define in the asymmetric system 
than in the symmetric one of In the symmetric system the unbound 



13 





-10 -5 




Figure 6: Plots of the trajectories for the upper skyrmions of a system with 
an asymmetric potential barrier of width 6=2 for various values of T and 
d : a) T = 0.3 and d = 4, b) T = 0.5 and d = 4, c) T = 0.5 and d = 5 and 
d) T = 0.8 and d — 5, where the horizontal axis represents the x coordinate 
and the vertical axis the y coordinate. 



14 



skyrmions, after being repelled from the barrier, were unable to move into 
the lower or upper semi-planes because the obstruction was placed over the 
whole x range. In the asymmetric system the skyrmions are able to move 
into these semi-planes. Thus in the symmetric system the trajectory of an 
unbound configuration was clearly defined. When two skyrmions are close 
enough to interact they will perform their usual circular motion about the 
configuration centre. It is difficult therefore to distinguish between states 
that are unbound or those that are not. The only way one is certain the 
skyrmions are unbound is to observe states for a sufficiently large values of T. 
In those cases, after the skyrmions have moved around the the obstructions, 
we observe the skyrmions moving about the boundary of the system singly. 
Periodically the skyrmion configuration undergoes circular motion when the 
skyrmions are close enough to interact. 

Ill Potential hole 

In [5] it was observed that the dynamics of the skyrmion configuration in- 
teracting with a potential hole was largely dependent upon whether the 
skyrmions were bound or not. When the skyrmions were unbound they 
were free to separate from each other. During the separation the skyrmions 
moved along the edge of the potential hole. This movement continued until 
the skyrmions reached the physical boundary of the system. It is here that 
both skyrmions underwent reflection at right angles to the boundary, travers- 
ing the hole and then executing an identical trajectory along the opposite 
edge of the the hole. The skyrmions motion continued back to the starting 
points where the system looked the same, with just an exchange of the two 
skyrmions. In the system of an asymmetric obstruction we would expect the 
previously mentioned symmetric trajectories to no longer persist. As indi- 
cated in previous sections, the hole exists on only half of the plane and hence 
the dynamics in this system should differ from the symmetric system. How- 
ever, we expect that the factor determining the form of the dynamics should 
be the same as in the symmetric system, namely, the binding energy, Eb, 
defined in ( 1131) . In the following three sections we shall present and discuss 
the results of a two-skyrmion configuration interacting with a potential hole 
of varying widths and depths, for a number of different skyrmion separations, 
d. The different sections are categorised according to whether the skyrmion 
dynamics is bound, unbound or in between. 
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Table 2: Table showing the variation in the total energy and binding energy 
for a d = 4 two-skyrmion configuration interacting with a potential hole of 
width b = 2 for various values of T. 

IIL1 Bound dynamics 

The skyrmions are initially set up in an attractive channel, which means 
that for close distances they form a bound state. Thus the states that one 
would presume to exhibit bound-like behaviour are those corresponding to 
small distances of separation and small hole depths. This natural observation 
turns out to be correct and the skyrmions, for small values of the hole depth 
T and distance of skyrmion separation d, do behave as a bound state. Tables 
[2] and [3] demonstrate the variation in the binding energy for the skyrmion 
configuration with the hole depth, for a fixed value of the width b = 2. The 
data shows that as the hole depth gets larger, the binding energy increases. 
As we will discuss in the later sections, this method of analysis does not 
include the effect of the hole on the skyrmion during its motion. Moreover, 
the binding energy is more of a dynamic quantity. Skyrmions are extended 
objects and in an asymmetric system one of the skyrmions feels the effect 
of the hole more than the other. Although the energy of each skyrmion can 
vary, for small values of T and d the static binding energy provides a good 
explanation of the behaviour of the skyrmions. 

Fig. [7] shows the trajectories of the upper skyrmion of a d — 4 configura- 
tion interacting with a potential hole for various values of the hole depth T. 
The trajectories show that the skyrmions are still bound, as was suggested 
by the binding energy analysis. Each plot was taken for a time length of 
1000 sees. The behaviour of the skyrmions in this system is reminiscent of a 
sling-shot type motion. Due to the small value of d, the skyrmions are tightly 
bound. When the skyrmions initially try to execute the normal circular mo- 
tion, the upper skyrmion, s 1; is moving towards the hole and the lower, s 2 , 
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Figure 7: Plots of the trajectories for the upper skyrmions in a system with 
an asymmetric potential hole of width b = 2 with d = 4 for various values 
of T: a) T = -0.1, b) V = -0.2, c) T = -0.4 and d) V = -0.5, where 
the horizontal axis represents the x coordinate and the vertical axis the y 
coordinate. 
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is moving into the space away from the hole. As si approaches the hole its 
energy is reduced due to its tail interactions with the hole. To conserve the 
total energy of the system the skyrmions separate a little from each other. 
This separation increases as s\ gets into the hole. In the meantime S2 moves 
further away from the hole as s\ approaches it. Classically we know that 
objects in a conserved system should speed up as they traverse a potential 
hole. Thus as Si enters the hole it quickly traverses it. The interaction be- 
tween the two skyrmions during this period results in S2 being sling shot 
around the hole towards the boundary at infinity. Once both skyrmions are 
out of the hole the normal motion can resume and the distance between the 
skyrmions corrects itself, resulting in S2 being kept to the left of the hole, as 
we observe it. This motion continues as the skyrmions move along the edge 
of the hole. The plots corresponding to larger values of |T| show that the 
skyrmions have travelled further than for the smaller values of |T|, over the 
same time length. This is due to the hole depth, where for a larger value of 
|T| the skyrmions move faster across the hole, resulting in the increase in the 
distance travelled in the same amount of time. 

Although not shown here, the skyrmion configuration for d — 5 produces 
similar plots to Fig. [7J with a few differences. For a d = 5 system the 
skyrmions, although still executing circular motion about the configuration 
centre, follow a trajectory that is more distorted than the skyrmions of Fig. 
[7J The skyrmions in this system are also able to penetrate the hole more than 
the skyrmions in the d = 4 case. The ability to penetrate the hole is due to 
the skyrmions' binding energy. As before, let si be the upper skyrmion and 
S2 the lower. The skyrmions try to execute their normal circular motion but 
as s\ approaches the hole its energy decreases. The initial binding energies 
for the d = 4 and d = 5 system are largely different due to the attractive 
channel in which they were created. As the energy of one skyrmion decreases 
the skyrmions separate from each other to conserve energy. Consequently, 
the binding energy or the interaction between the skyrmions, is reduced. In 
the system for d = 4, the skyrmions do not separate as much from each other 
as in the d = 5 case due to the tightness of the bound state. In the d = 5 state 
the initial maximum separation between the skyrmions occurs at the point 
when si is in the hole. The weak interaction between the skyrmions slows 
down the distorted rotation about the configuration centre. In comparison 
with the trajectories of Fig. [7J the skyrmions' rate of rotation is not as 
reduced by the hole as in the d = 5 system. Since the skyrmions are much 
more tightly bound at the start, the effect of the hole upon the interaction 
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Table 3: Table showing the variation in the total energy and binding energy 
for a d = 5.5 two-skyrmion configuration interacting with a potential hole of 
width b = 2 for various values of T. 



between the skyrmions is much less pronounced. Therefore the skyrmion 
trajectory is less distorted and the skyrmions move faster while traversing 
the hole. 

III. 2 Transition dynamics 

In this section we describe the transition from the bound systems described 
in the previous section to one in which the skyrmions behave as if they were 
unbound, similar to the dynamics seen in the symmetric system of [5j. The 
transition state is best observed through the variation in T for the d = 5.5 
system. It can be seen from the trajectories in Fig. [8] that the skyrmions for 
T = —0.2 still form a bound state. However, as we increase the hole depth 
from T = —0.2 to T = —0.3 and further to T = —0.5, the dynamics is no 
longer bound and the skyrmion trajectories are completely uncorrelated, as 
seen in Fig. [91 In Fig. [H^i and b the skyrmions are able to penetrate the hole 
more than in other previous systems. As one skyrmion penetrates the hole 
the other remains stationary outside. The skyrmions in this state are just on 
the point of being broken up by the hole. The skyrmions still remain bound 
after the upper skyrmion comes out of the hole and this is evident in the 
trajectory. In the case of T = —0.4, —0.5, the skyrmions behave differently. 
Here one of the skyrmions moves along the axis of the hole after a small initial 
circular motion. The other comes to rest after this initial movement. The 
other skyrmion continues along its path until it reaches the boundary, where 
it gets reflected and moves back towards the stationary skyrmion. When 
they get close to each other, the skyrmions violently oscillate due to their 
mutual interaction and these oscillations are discussed in a later section. 
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Figure 8: Plots of the trajectories for the upper and lower skyrmions of a 
system with an asymmetric potential hole of width 6 = 2 with d = 5.5 for 
various values of T: a) T = —0.3, upper skyrmion , b) V — —0.3, lower 
skyrmion, c) T = —0.2, upper skyrmion and d) F — —0.2, lower skyrmion, 
where the horizontal axis represents the x coordinate and the vertical axis 
the y coordinate. 
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Figure 9: Plots of the trajectories for the upper and lower skyrmions of a 
system with an asymmetric potential hole of width 6 = 2 with d = 5.5 for 
various values of T: a) T = —0.5, upper skyrmion, b) T — —0.5, lower 
skyrmion, c) T — —0.4, upper skyrmion and d) F — —0.4, lower skyrmion, 
where the horizontal axis represents the x coordinate and the vertical axis 
the y coordinate. 
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Let us consider the binding energy, Eb, defined in (1131) to be a function 
of D, the distance a single skyrmion is from the asymmetric hole in (1131) . By 
examining the binding energy for the skyrmions in these systems as a function 
of D, it can be seen how the skyrmions become unbound as they get close 
to the hole. By examining tables [2] and [3] and utilising other information, 
we can examine how the single skyrmion energy varies with D. The energies 
of the two skyrmion configuration in the presence of a hole of depth T = 
-0.4 is E 2 = 2.0977/8tt and for T = -0.5 is E 2 = 2.0936/8tt. Tables 1 
and [5] show how the binding energy of the skyrmion configuration varies 
when one of the skyrmions approaches the hole. It is clear that when the 
skyrmions, in the case of T = —0.4, get near D = 1 the binding energy of 
their configuration becomes positive and the skyrmions are no longer bound. 
Thus as the skyrmion approaches the hole, after the initial circular motion, 
its starts to feel the effect of the hole. To conserve energy the skyrmions must 
increase their distance of separation. In order to do this the upper skyrmion 
moves along the edge of the hole. The more the skyrmion feels the hole 
the more the configuration separates. As the binding energy increases the 
skyrmions interact less and the upper skyrmion is free to move along the edge 
of the hole towards the boundary. After the reflection from the boundary the 
skyrmion speeds up traversing the obstruction, moving towards the opposite 
edge of the hole. This process then repeats itself; when the skyrmion reaches 
the stationary one it interacts with it causing the the stationary skyrmion to 
begin moving. The previously moving skyrmion now comes to rest. It was 
evident in the previous section that as T increases, for a given value of d, the 
skyrmions' trajectories change. The skyrmions cease to behave as a bound 
state and the behaviour of the two skyrmions becomes uncorrelated. 

The initial distance of closest approach of the upper skyrmion can be 
estimated from our dynamic binding energy arguments. In the T = —0.4 
system, the closest distance the skyrmion reaches, based upon the trajectory, 
is approximately D = 1.3. Comparing this with table IH we can see that 
the transition point occurs between D = 1.5 and D = 1. Based upon the 
numerical values of each, the transition point is closer to 1.5 than 1. In 
the T = —0.5 case, the closest distance is about D = 1.75, based upon the 
trajectory. Estimating this distance from the energetics, we see that the 
transition point occurs between D = 2 and D = 1.5, with it being nearer to 
1.5. Thus from just energetic arguments we have a reasonable estimate of 
the closest distance the skyrmion can initially get to the hole. 

In the section on the potential barrier system we examined the effect 
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D 


£i/8tt 


Eb/Stt 


2.5 


1.0623 


-0.0269 


2.0 


1.0566 


-0.0155 


1.5 


1.0503 


-0.0029 


1.0 


1.0458 


+0.0061 


0.5 


1.0432 


+0.0072 



Table 4: Table showing the variation in the binding energy for a T = —0.4 
hole of width 6 = 2 for various values of the single skyrmion distance from 
the hole D. 



D 


£i/8tt 




2.5 


1.0613 


-0.029 


2.0 


1.0534 


-0.0132 


1.5 


1.0455 


+0.0026 



Table 5: Table showing the variation in the binding energy for a T = —0.5 
hole of width b = 2 for various values of the single skyrmion distance from 
the hole D. 
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of increasing the bounding box of the simulation. It was found that the 
skyrmions trajectory away from the obstruction was uniquely determined 
by the boundary. In the potential hole system the boundary once again 
determines the dynamics of the skyrmions. The skyrmions in such a system 
often reach the boundary where they get reflected. If identical skyrmions are 
put into a system with a larger grid, the skyrmions undergo a reflection from 
the boundary at a later time than for a system of a smaller grid. 

IV Angular Momentum 

In this current study we have considered the scattering of a skyrmion config- 
uration off an asymmetric obstruction, similar to the symmetric obstructions 
of [5], except now the obstruction no longer exists for x > 0. In this sec- 
tion we shall derive the contributions the potential obstructions make to the 
derivative of the orbital angular momentum I. We shall also discuss the 
trends seen in the data obtained for the orbital angular momentum, /. It will 
be seen that the time variation in I can be decomposed into the behaviour 
of the average size of the skyrmions, r, and the modulus of the skyrmion 
guiding centre, R, both defined previously. It was noted in the potential 
hole section that the simulations of the skyrmions in the transition dynam- 
ics exhibit oscillations in the energy density. Here we shall show that these 
oscillations play a significant role in the conservation of the orbital angular 
momentum. 

IV. 1 Asymmetric potential obstructions contribution 
to 1 

The potential obstructions are introduced as an inhomogeneity in the po- 
tential term's, V(4>), coefficient 73. The obstructions can be defined more 
concisely if they are written using Heaviside functions. In this formalism we 
can rewrite the approximation of the obstruction as: 

v{±) = i(i - 4>l) + ^r(i - 4){e(-x + x Q ) [e(y + yo ) - e(y - y )] }. 
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We now proceed, as in to calculate the contribution the obstruction makes 
to q: 



lobs 



-e„ud, 



3 ^0 3 (1 + Te(-x + x ) Q(y + y )-Q(y-y ) )}. 



The time derivative of the orbital angular momentum was defined previously 
by ffTUl) . If the above is inserted into ffTUj) . then an expression for the contri- 
bution that the asymmetric potential obstruction makes to I can be found: 



obs 



j2 2 ■ 

d xx q obs 



d 2 xx 2 d y (-<J)l) \5(-x + x ) [&(y + y ) - &(y - yo) 
d 2 xx 2 d x (^4>l)Q(-x + x )5(y + y ) 



r f i 

- / d 2 xx 2 d x (-(j) 2 3 )e(-x + x )5(y-y ). 



These integrals can be easily computed to give the overall contribution of 
the asymmetric obstruction to I, given by: 
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dy xjd,, 
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dy x 2 d, 
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yy^' 



X=XQ 



j2 2 • 

d xx q ohs . 



(14) 



IV. 2 Examination of J and J 

In any system the total angular momentum, J = I + m, should be con- 
served. However, it was discovered in [5] that for scattering processes in 
Landau-Lifshitz models the total angular momentum, as defined by ([5]), is 
not conserved. On further examination it was found that throughout all the 



25 



simulations, the total magnetization in the third direction, m, is well con- 
served and it is the orbital angular momentum, I, which is not conserved with 
time. In the asymmetric system such a situation presents itself again, where 
the total angular momentum is not conserved due to the non-conservation 
of the orbital angular momentum (jSJ), denoted I fields- The non-conservation 
of J is apparent in the plots of Fig. [TUJ Both of these plots show the total 
angular momentum, J = I + m and its components, for a d = 6 skyrmion 
configuration, interacting with a potential hole or barrier with T = ±0.25. 
Here it is evident that for either a hole or a barrier the total angular mo- 
mentum is not conserved due to the non-conservation of /, since m is well 
conserved throughout. Before analysing this apparent non-conservation of J 
we shall firstly discuss the general behaviour of / and hence J, during the 
simulations. 

The time evolution of I, as formulated by Papanicolaou and Tomaras [7J, 
provides us with an insight into the behaviour of the skyrmions during the 
simulations. It was already seen that the average skyrmion size, r and the 
guiding centre coordinate, R, can be related to the orbital angular momen- 
tum, /, by : 



or 

/ = 2nQ{r 2 + R 2 }. 

The guiding centre coordinate, R, for a single skyrmion system provides 
the location of the skyrmion centre. In the case of a two-skyrmion configu- 
ration, R provides the position of the centre of the configuration. The initial 
placement of the skyrmions results in the centre of the configuration nearly 
positioned at (0,0). In the case of a hole or barrier, the asymmetry of the 
system results in the position of the configuration centre becoming a function 
of time. How it varies during the simulation depends on the system. In the 
case of a potential barrier the skyrmions collectively move about the system. 
Therefore as the skyrmions begin to move, R deviates from (0, 0) and \R\ 
steadily increases. \R\ continues to increase until the skyrmions reach the 
boundary of the system. As the skyrmions orbit the boundary, \R\ slowly 
increases and decreases as the skyrmions move towards and away from the 
boundary. If we consider the behaviour of I as corresponding to the dual 
behaviour of the quantities r and R, then as \R\ steadily increases, so does 
I. This overall trend in \R\ is evident in Fig. [TOa which clearly shows / 
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Figure 10: Plots of the total angular momentum, J = I + m, and its com- 
ponents as a function of time for a d = 6 skyrmion configuration interacting 
with a potential hole or barrier for T = ±0.25: a) hole and b) barrier. 
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steadily increasing with time. In the case of a potential hole we can use 
similar analysis. Fig. [TUb shows a plot of the orbital angular momentum for 
potential hole with T = 0.25. Here the increase in Z, is much larger than in 
the case of the barrier system. In this particular system, the configuration 
centre moves much more than the barrier system and is periodic. This is 
evident if we recall the dynamics of the skyrmions in this system. Here one 
of the skyrmions moves along the edge of the hole while the other remains 
stationary. Thus the centre of the configuration moves away from (0, 0) as 
the skyrmion approaches the boundary. After the skyrmion is reflected, R 
begins to approach (0, 0) once again. However, due to the dynamics R does 
not reach (0, 0) again. This periodic increase and decrease in \R\ matches 
the behaviour of I for the potential hole. 

The behaviour of the orbital angular momentum, with R — 0, can be 
understood exactly in terms of the behaviour of r. The skyrmions tail is 
governed by value of the potential coefficient, 73 and the regions where this 
is modified means the skyrmion tails become modified. As the skyrmion 
moves into a region for which 73 — > 1 + T, its tail contracts or expands. This 
continued shrinking and expansion, as the skyrmion moves on and off the 
potential obstructions, results in a time oscillation of r 2 . In Fig. [TUh the 
behaviour of r is evident in the time dependence of I. The periodic small 
oscillations occurring in I, correspond to the time variation of r. It was 
seen in the symmetric barrier system that as the skyrmion configuration is 
'on' the obstruction, the skyrmions reduce their distance of separation. In 
the asymmetric barrier system the small oscillations in I are caused by this 
behaviour. The interplay between the time oscillation in r 2 and the gradual 
increase in R? sufficiently accounts for the observed behaviour of the angular 
momentum. 

It has already been noted that the total angular momentum, J, is not 
conserved. It was discovered in our previous work that by examining how 
the potential obstructions were constructed, an analytical expression for the 
contribution the potential obstructions make to the time derivative of the 
orbital angular momentum, /, was found. These expressions were in the 
form of two integrals which were calculated during the simulations. When 
the contribution of the potential obstructions was included in the expression 
for J, the conservation of total angular momentum was restored. These 
same integrals, for the asymmetric obstruction, were computed in section 
IIV.1I and were calculated throughout the simulations. The contribution the 
obstruction makes to the derivative of the total angular momentum, J tot , is 
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denoted l obs . The contribution to J to t from the fields, given by is denoted 
I fields- Therefore the overall conservation of the total angular momentum can 
be written: 

Jtot lobs I fields i 

where there is no contribution due to the total magnetization in the third 
direction, m, since throughout rh = 0. Fig. [TT] shows the plots of the total 
angular momentum, J tot = l bs + 1 fields for two different systems. In Fig. [TTh . 
lobs, I fields and Jtot are plotted for a d — 6 skyrmion configuration interacting 
with a potential barrier of width 6 = 2 and height T = 0.25. It can seen 
from the plot that the overall behaviour of I fields and l b s match, resulting 
in J tot being well conserved in time. It can be seen in the latter stages of 
the plot that this conservation is slightly destroyed. The skyrmions at this 
point start to appreciably interact with the boundary of the system. It is 
due to these boundary effects that J to t ^ at these points. Fig. [TTb shows 
the corresponding plot for a d = 5 skyrmion configuration interacting with 
potential hole of depth T = —0.5. It is evident here that the behaviour of 
I fields and l bs is again nearly identical, resulting in the overall conservation 
of Jtot- We note again that there exist regions in Fig. [TTb where J to t is not 
conserved. In this particular system these regions correspond to the points 
when one of the skyrmions is reflected from the boundary and constitute 
another boundary effect. 

IV. 3 Skyrmion oscillations 

It was seen in the previous section that the total angular momentum, J to t, 
for the asymmetric system is well conserved, just as it was for the sym- 
metric system. There were regions where this conservation was temporarily 
destroyed and we have noted that such effects arise due to the skyrmions 
interaction with the boundary of the system. However for critical values of 
T and d, during the interaction of a skyrmion configuration with a potential 
hole, there exists regions where Jtot is not conserved. Here it shall be seen 
that the regions where J to t 7^ do not arise due to to boundary effects but 
are due to an important physical process. 

Fig. [T2] shows the plot of J to t for a d = 5.5 skyrmion configuration 
interacting with potential hole of width b = 2 and depth T = —0.5. The 
trajectories of both the upper and lower skyrmions in this system are plotted 
in Fig. [13]over the period of time t = 500—1000 sees. This period corresponds 
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Figure 11: Plots of the time derivative of the total angular momentum, J to t 
and its components I fields and l f, s as a function of time for a two different 
systems: a) d = 6 skyrmion configuration with b = 2 for T = 0.25 and b) 
d = 5 skyrmion configuration with b = 2 for T = —0.5 
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to the same period in Fig. Q2J when the large spike in J tot is observed. The 



Figure 12: Plot of the derivative of the total angular momentum, J to t = 
I fields + Lbs for a d = 5.5 two-skyrmion configuration interacting with a 
potential hole of width 6 = 2 and T = —0.5. 



points P,Q and R of Fig. [T5] can be found by examining Fig. which 
show plots of y(t) for both the upper and lower skyrmions. The locations of 
the upper and lower skyrmions for the points P,Q and R are shown on the 
trajectory in Fig. [13] and are given the subscript 1 and 2 to identify the upper 
and lower skyrmions respectively. During the periods of P-Q and Q-R, Jtot 
undergoes a rapid change. Over the same time length both the skyrmions 
have traversed the hole, firstly the upper skyrmion, si, and then the lower, 
S2- In the path, Pl-Ql, si traverses the hole while s 2 starts to move along 
the edge of the hole. In the next path, Ql-Rl, si now moves along the lower 
edge of the hole and begins to execute circular motion after its interaction 
with s 2 . The path Q2-R2 coincides with s 2 traversing the hole. We must 
therefore ask ourselves, what distinguishing features are present during the 
time periods P-Q and Q-R that are absent the rest of the time? To do this we 
examined a simulation of the energy densities of the skyrmion configuration 
for the d = 5.5, T = —0.5 system, during this time period. Fig. [15] shows 
snap shots of the energy densities of the upper and lowers skyrmion in this 
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system. The snap shots evolve from a) to d) in 10 second increments. In such 
a system one can see the skyrmions undergo large oscillations in their energy 
density during their motion. In the simulation of this system, s\ traverses 
the hole during the initial 500 sees. After s\ traverses the hole for a second 
time, the skyrmions undergo a violent oscillation in their energy densities. 
Their energy densities dramatically pulse before the skyrmions continue with 
their motion around the edge of the hole. The first of these oscillations occur 
during the period P-Q and the skyrmions are subsequently seen to undergo 
another oscillation, during the period Q-R. The large change in the energy 
densities, causes the observed spike in Fig. [121 There exist other examples of 
skyrmions oscillating in these systems and in all cases the conservation of J to t 
is destroyed during such periods. This evidence leads us to the observation 
that the non- conservation of Jtot and skyrmion oscillations in Landau-Lifshitz 
models are intertwined. 



V Conclusion 

In the asymmetric barrier system the skyrmions behave similarly to the 
skyrmions observed in the symmetric barrier system. The bound skyrmions 
perform the normal circular motion about the configuration centre moving 
on and off the barrier during their evolution. The distance of separation of 
the skyrmions increases and decreases with time. The maximum separation 
between the two skyrmions corresponds to the times when one or both of the 
skyrmions is off the obstruction. The minimum separation corresponds to 
times when both the skyrmions are interacting appreciably with the barrier. 
This change in the distance of separation between the skyrmions results in 
the distortion of their path as they traverse the obstruction. In this same 
system the skyrmions are able to move into the upper and lower planes. In 
the symmetric system they could not do this due to the obstruction existing 
over the whole x range. When the skyrmions have moved away from the 
obstruction they orbit the boundary of the system. It was noted that this 
phase of the dynamics is uniquely determined by the boundary of the sys- 
tem. If the skyrmions are placed in an identical system, with a larger grid, 
after initially executing identical trajectories as the skyrmions for the system 
of a smaller grid, the skyrmions move in a near identical motion about the 
boundary of the system, except they are further from the obstruction than 
the skyrmions of a smaller grid. The binding energy in the barrier system 
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Figure 13: Trajectory of both skyrmions for ad— 5.5 configuration inter- 
acting with a potential hole of width b = 2 and depth T = —0.5, plotted 
over the time range 500 — 1000, where the horizontal axis represents the x 
coordinate and the vertical axis the y coordinate. 
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Figure 14: Plots of y(t) for the upper (b) and lower (a) skyrmions in a system 
with an asymmetric potential hole of width 6 = 2 and T = —0.5. 
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Figure 15: Plots of the energy densities of the upper and lower skyrmions for 
a d — 5.5 skyrmion configuration interacting with a potential hole of width 
6 = 2 and depth T = —0.5. The time line for this set of plots begins with a) 
through to d) in 10 second increments. 
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also plays an important role. The bound state formed by two skyrmions can 
be broken by a potential barrier. The breaking point in the asymmetric sys- 
tem is much more difficult to observe than in the symmetric system. In the 
symmetric system the geometry of the potential obstructions uniquely deter- 
mines the trajectories of the unbound skyrmions. The unbound skyrmions in 
the asymmetric system move almost independently of each other about the 
system. Periodically the skyrmions, when they are close enough to interact 
due to the confinement of the boundary, execute circular motion about their 
centre. 

In the system of a potential hole the skyrmions behave completely dif- 
ferently from the skyrmions observed in the symmetric system. Here the 
binding energy plays a crucial role in categorising the different phases of the 
dynamics; bound, unbound or in between the two. The skyrmions, for tightly 
bound configurations, orbit each other as they traverse the edge of the hole. 
These bound states occur for small distances of separation and hole depth. 
As the hole depth or the distance of separation increases, the binding energy 
of the skyrmions increases. However, as the skyrmions approach the point at 
which the bound state is broken, the dynamics of the skyrmion configuration 
changes. The skyrmions in these 'in between states' are able to penetrate 
the hole. They behave as bound skyrmions in the sense that their trajectory 
could be described as a very distorted circular path but they behave as un- 
bound skyrmions since they separate appreciably during the simulation. In a 
potential hole system, as was noted in the analysis on the symmetric system, 
the binding energy becomes a dynamic quantity. As the skyrmions approach 
the hole the binding energy or the interaction energy between the skyrmions, 
changes. The skyrmions are extended objects. When one of the skyrmions 
approaches the hole, the energy of the skyrmion decreases due to the tail 
interactions with the hole. Due to the asymmetry of the system, one of the 
skyrmions can feel the effect of the hole more than the other. Energy conser- 
vation allows the skyrmions to separate. The more a skyrmion approaches 
the hole the more they can separate and thus the skyrmion bound state is 
broken by the potential hole. This process is sensitive to the hole depth, T. 
In the transition state dynamics, d = 5.5, increasing the value of T results 
in a complete change in the behaviour of the skyrmions. This is largely due 
to the binding energy as the skyrmion initially approaches the hole. In the 
case of large values of T, the skyrmions know they are unbound when they 
are further from the hole and begin to execute dynamics reminiscent of the 
skyrmions for a symmetric potential hole. When V is slightly smaller, due to 
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Eb-i the skyrmion is able to penetrate the hole and there is still an apprecia- 
ble interaction with the other skyrmion. This allows the skyrmion to escape 
the hole and a pseudo-bound state behaviour to resume. Our examination of 
the dynamic behaviour of the binding energy led us to a reasonably accurate 
estimate for the closest initial distance of approach of the upper skyrmion 
during the transition dynamics. We thus believe that the subsequent analysis 
regarding the binding energy and the skyrmions' dynamics dependence upon 
it, is a valid and a good representation of the details of these complicated 
systems. 

The system of unbound skyrmions in the potential hole system also shows 
some unusual properties, some of which were seen in the transition dynamics. 
The unbound skyrmions occur for large values of the distance of skyrmion 
separation and large hole depths. The skyrmions of these system behave, 
as was alluded to previously, similarly to the skyrmions in the symmetric 
potential hole system where both skyrmions moved along the axis of the hole. 
Here, due to geometry of the hole, only one of the skyrmions moves along 
the edge of the hole, while the other, after the initial motion, comes to rest. 
During the initial period of separation both of the skyrmions' energy densities 
are seen to oscillate. This oscillation continues until the skyrmions reach a 
critical distance of separation, where it ceases. These oscillations surface 
again after the skyrmion is reflected from the boundary, when it has reached 
another critical distance from the stationary skyrmion. These oscillations 
are rather peculiar and the only other system to bear such oscillations is the 
transition dynamics. In the transition dynamics for V = —0.4 and —0.5, the 
skyrmions' energy densities undergo violent oscillations. These oscillations 
are the result of the excitation of one or more of the internal modes of the 
skyrmions due to the energetics during the skyrmion separation. We would 
expect that any of the excited modes would be of the lowest order. Piette 
et al [8] have previously undertaken an analysis of the excited modes of the 
skyrmions of the new baby skyrme model, in the fully relativistic system. In 
[8] the modes of a single skyrmion and a two-skyrmion configuration during 
a scattering process were analysed. It would be interesting to see how the 
excited modes of the skyrmions observed in our simulations compare with 
the excited modes seen in their work. 

One of the most interesting observations of our simulations was the ap- 
parent non- conservation of the total angular momentum J (given its usual 
definition). This non-conservation of J was due to the non-conservation of 
the orbital angular momentum, I, as we have found that in all of the simu- 
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lations the total magnetization in the third direction, m, was well conserved 
in time. At the same time we showed that I ^ 0. Thinking about this fur- 
ther we showed that when a system involves potential obstructions, these 
obstructions made a significant contribution to /. Hence one has to modify 
the conventional definition of I. We have found and calculated this missing 
contribution. In all of the simulations of the skyrmions in both the symmetric 
and asymmetric systems we have shown that the addition of the contribution 
from the obstructions compensates I. This results in the overall conservation 
of I and J for the full system. However, there were a few systems studied 
that didn't adhere to the above statement. In the transition dynamics in the 
asymmetric system, there were regions where the angular momentum due 
to the fields, I fields, and the angular momentum due to the obstruction, l obs , 
did not sufficiently cancel each other out, in order for total angular momen- 
tum conservation. However, when comparing the size for which J 7^ in 
such regions with the total angular momentum as a whole, these regions of 
non-conservation are brief and the magnitude in most of the cases is small. 
In some systems, where the non-conservation of J was larger, we noted that 
this was due to boundary effects. 

An examination of the times when there was an appreciable discrepancy 
in I fields + lobs 7^ 0, provided a direct link between the skyrmion oscillations 
and the non-conservation of J. Initially J is well conserved but as the os- 
cillations begin to increase this conservation is destroyed. There reaches a 
point when the J conservation is restored. It is at these points that the 
skyrmion oscillations have ceased. In many of the cases these oscillations 
reappear again and again during the skyrmions dynamics and each time the 
conservation of J is destroyed. The size of the oscillation and the magnitude 
to which J ^ are also correlated. In the transition dynamics the large 
spike observed in the T = —0.5 plot of I fields + lobs system corresponds to 
the violent oscillation during the dynamics of the skyrmions. Conversely in 
the system with T = —0.3 smaller oscillation are apparent in that system, 
resulting in a smaller spike in I fields + lobs- The regions where J to t ^ sug- 
gest that the inclusion of an additional term is required for the overall exact 
conservation of the total angular momentum, J, for skyrmion scattering in 
a Landau-Lifshitz model. Any term or terms that are included must have 
a direct correspondence with the oscillations of the skyrmions observed in 
these systems. 
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VI Appendix: Numerical procedures 



Unfortunately, it is impossible to solve ([3]) analytically we have therefore 
had to study this problem numerically. The fields and their derivatives were 
discretised in the usual manner and were placed on a lattice of 251 x 251 
points, with lattice spacing dx = 0.1. The numerical integration of the 
3-coupled differential equations of ([3]) involved the use of a 4 th order Runge- 
Kutta method of simulating time evolution with a time step of dt = 0.001. 
The various integrals calculated throughout the simulations were performed 
using a 2-D Simpson's rule. The constraint equation requires that the fields 
lie on the 2-sphere, (j) =1, and this was imposed at every time step by 
rescaling each field component so that fa — > fi— . 



The skyrmions were initially placed at (0, ±d/2) in the upper and lower 
planes, where d is the distance between the two skyrmion centres (xi,yi). 
The trajectory of each skyrmion was tracked by following the maxima of 
the topological charge density and interpolating between the lattice points. 
The coefficients 7; have been set to unity in all the simulations unless stated 
otherwise. 
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